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By means of a projector-operator formalism we derive an approximation based on a self consistent hybridization 
expansion to study the ground state properties of the Anderson Impurity model. We applied the approximation to 
the general case of finite Coulomb repulsion U, extending previous work with the same formalism in the inhnite-[/ 
case. The treatment provides a very accurate calculation of the ground state energy and their related zero temperature 
properties in the case in which U is large enough, but still finite, as compared with the rest of energy scales involved 
in the model. The results for the valence of the impurity are compared with exact results that we obtain from equations 
derived using the Bethe ansatz and with a perturbative approach. The magnetization and magnetic susceptibility is also 
compared with Bethe ansatz results. In order to do this comparison, we also show how to regularize the Bethe ansatz 
integral equations necessary to calculate the impurity valence, for arbitrary values of the parameters. 
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1 Introduction The Anderson impurity model (AIM) 

ffl , which is one of the most studied Hamiltonians includ¬ 
ing strong correlations in condensed-matter physics, has 
been solved exactly by means of the Bethe ansatz |2][3]]. 
In addition, the spectral function and many other correla¬ 
tion functions have been accurately computed by using the 
numerical renormalization group (NRG) II41I51. 

However, several approximated schemes are also used. 
Frequently, they permit a rapid scanning of the properties 

as a function of the model parameters, they shed light to the 
expected behavior and can be used to fit experiments, when 
due to a large number of degrees of freedom it is not possi¬ 
ble to use more robust but also time consuming techniques. 
Moreover, there is no Bethe ansatz solution in the case of 
a frequency dependent hybridization function. Among the 
approximations, the self consistent hybridization (SCH) 
expansions for solving the AIM are commonly used within 
a large class of different problems. The non-crossing ap¬ 
proximation, NCA, j6j which represents the simplest fam¬ 


ily of these self consistent treatments, provides an accurate 
calculation of the Green functions, as well as many other 
properties, when the Coulomb repulsion is taken to be in¬ 
finite. Since the 80’s the NCA has been successfully ap¬ 
plied to study Ce compounds 0, non-equilibrium trans¬ 
port properties ED and, more recently, non Fermi liquid 
behaviors 0, electron-phonon interaction ca and spec¬ 
troscopy of a double quantum dot system DU- 

When the Coulomb repulsion U takes a finite value, 
the NCA has failed to give the correct Kondo scale ( Tk )■ 
Unfortunately its value is found to be very underestimated 
as compared with the correct one obtained from the Bethe 
ansatz solution of the model. The next leading order in the 
self consistent expansions that partially solves this pathol¬ 
ogy is often known as the one crossing approximation, 
OCA, HT2HT31H41 . Within this extended formalism other 
classes of problems have been successfully investigated 
E). Among them, its mayor application is in the context 


Copyright line will be provided by the publisher 







2 


I. J. Hamad et al.: Ground State of the AIM 


of the dynamical mean-field theory as an impurity solver 

ED- 

However, the use of NCA or OCA approximations im¬ 
poses some limitations. Finite temperatures T have to be 
used, specifically temperatures larger than 0.01 Tk in the 
former and 0 . 114 : in the later, in order to avoid an arti¬ 
ficial increase by abut 10% in the spectral weight at the 
Fermi level. The same pathology arises when these tech¬ 
niques are applied to systems in which the ground state 
without hybridization is a non degenerate state. This is the 
case when either the impurity state is empty (or formed by 
an even number of electrons) or when a magnetic field is 
applied to the impurity that brakes the Kramers degeneracy 
(Zeeman effect). 

Remarkably, these limitations when calculating dy¬ 
namical properties, are absent in the case of static ones IfTTl 
like those derived from the ground state energy. In the 
early approach given by InagakiQjO and Keiter et «/. ||19l 
the valence of the impurity as well as the charge and spin 
susceptibilities were calculated from the ground state en¬ 
ergy (T = 0) in the U infinite limit and for a large value of 
the conduction bandwidth D —> oo. 

Recently, an improvement of these calculations incor¬ 
porating higher order processes when calculating the oc¬ 
cupancy and magnetization of the impurity was done by 
three of us l20l . The agreement with exact results was re¬ 
markable. 

However, an extension of this static approximation to 
finite values of U is still desirable. This is justified due 
to the existence of relatively not so large values of the 
Coulomb interaction in bulk systems as well as in sys¬ 
tems of semiconductor quantum dots or carbon nanotubes 
through which transport experiments are currently being 
performed 12 1113111321 . 

In this contribution we extend the use of the SCH ex¬ 
pansions when calculating the ground state properties of 
the single impurity Anderson model to incorporate large 
but finite values of the Coulomb repulsion. We compare 
our approximated ground state properties, like valence and 
magnetization of the impurity with exact Bethe-ansatz cal¬ 
culations. We also provide a method to regularize the di¬ 
vergences that appear in some Bethe-ansatz expressions 
which might be useful for researchers interested in using 
these exact results. In addition, in the opposite limit, we 
start the discussion showing the excellent agreement of im¬ 
purity occupancy calculated from perturbation theory (PT) 
to second order in U, modified to give the correct result 
in the atomic limit and to satisfy the Friedel sum rule 1221 
|23| . for not so large values of U as compared with the hy¬ 
bridization strength. 

The paper is organized as follows. In Sec. [2] we in¬ 
troduce the model and, starting from a projector-operator 
formalism, we review the self consistent hybridization ex¬ 
pansion including finite values of the Coulomb repulsion, 
we generalize the self consistent hybridization expansion 
to include finite values of the Coulomb repulsion. In Sec. [3] 


we present and discuss the numerical evaluations of differ¬ 
ent physical magnitudes, including a detailed benchmark 
between Bethe ansatz results and the approximated PT and 
SCH-expansion. Finally, in Sec. [4] some conclusions are 
drawn. In the appendix we explicitly give the steps that we 
follow in order to regularize and compute in an efficiently 
numerical way the integral equations of the Bethe ansatz 
when calculating the impurity valence for arbitrary values 
of the model parameters. 

2 Model and formalism The Anderson Hamiltonian 
describes a system in which a localized interacting level is 
coupled with a conduction band of electrons, 

H = ^2 e k n ka + ^ E dndo + Un d ^n d |+ 

k * a (1) 

^(Vkdlcka + H.c.), 

ka 

where n ka = c k(T c kcr * s t ^ le number operator and c' ktT cre¬ 
ates a conduction electron with momentum k and spin cr 
in the conduction band. The operators n da = d\d,j and 
dj. refer to an electron in the interacting localized impurity 
state characterized by an energy E d and Coulomb repul¬ 
sion U. The coupling of this level with the bath is taken 
into account by the hybridization function given in terms 
of 14 by A(u) = n J2 k V£S{lo - e k ). 

We will be interested in looking for the energy of the 
ground state, Eo, of the system of A'-particles satisfying 
the eigenvalue equation, 

= E\i/j 0 )- (2) 

where |-0o) is the ground wave function of the system. 
From it, we can obtain other quantities such as the occu¬ 
pancy of the impurity, magnetization and susceptibility by 
differentiation with respect to E d or to the magnetic field 
B, which can be easily incorporated in this formalism via 
a Zeeman term gg,BB(n\. — n-p)/2 in the Hamiltonian. 

The procedure we follow in general is the same as 
the one applied in a previous work by three of us ll20l . 
The main idea is to divide the Hilbert space in two sub¬ 
spaces, by means of projector operators, but with the pe¬ 
culiarity that the working subspace consists of a single 
state. For that purpose we define the projector operator 
p i = n k *< kF n ka , where fci? is the Fermi wave vector 
and P 2 = 1 — P\. The action of the projectors is then 
to split the Hilbert space of A'-particles into two disjoint 
subspaces, namely S± which consists of only one state, 
|^i) = rifccr</c F C !crl0)’ the Fermi sea with no electrons in 
the (/-level (which is the ground state in absence of the im¬ 
purity) and S ‘2 with 1 02 ) representing the rest of the states 
in the Hilbert space. In principle, the state |<^i) can be cho¬ 
sen as the one in which the impurity is empty, singly or 
doubled occupied, arriving at the same final results. 


Copyright line will be provided by the publisher 





pss header will be provided by the publisher 


3 


The renormalized Hamiltonian that operates on the 
S i subspace can be obtained by eliminating |</> 2 ) of the 
Schrodinger equation (0, in a usual procedure [[26]. The 
action of the effective Hamiltonian on is given by, 

H\<h) = (Hu + H ^ e ^ Ho2 H ^) l^i) = E\<k), 

(3) 

where // tJ = T\HP 3 . The only contribution to the Hamil¬ 
tonian H 21 , that connects the subspace S\ with S 2 , comes 
from the hybridization term in Eq. 0. Its operation de¬ 
stroys an electron of the Fermi sea and promotes it into 
the impurity generating a singlet state belonging to the S 2 
subspace. 


principle, it is possible to obtain all the matrix elements, 
gij, of the resolvent operator G = (E — H 22)- 1 ■ However, 
to obtain the ground state energy, given in 0, it is nec¬ 
essary to calculate the diagonal and non diagonal matrix 
elements g k k' corresponding to the states below the Fermi 
level. The contribution to the energy of the non-diagonal 
path that starts from a hole in k and finishes with a hole 
in k! is of the order of 0(V 2 / 2D 2 ) times less than the di¬ 
agonal contribution. Following closely the treatment done 
in ]20l , and working in the thermodynamic limit, we begin 
by examining only the diagonal contribution g k k, which 
is expressed as a continuous fraction that ultimately, in the 
thermodynamic limit, can be written in a closed expression 
as: 


where the factor \/2 is included in order to normalize the 
new state. Here and in what follows, we employ the no¬ 
tation k and K for wave numbers below and above the 
Fermi level respectively. We emphasize that the state \4> k d) 
is, by construction, a singlet and therefore its total spin is 
St = 0. Applying (<j> 1 1 to the left of Eq. 0 we obtain an 
equation for the energy E, 


E = 6t + 2- 


N 


Y(m e n \<t>k'd), 


(5) 


k,k' 


where ex = 2 ]T ; , £fc represents the energy of the ground 
state of the Fermi sea. Also, in 0, we neglect the k de¬ 
pendence of the hopping 14 = V/\f~N as usual, where the 
scaling by 1 /\/N is necessary to obtain finite contributions 
to the energy of the Fermi sea-impurity hopping term. 

It can be noticed that to get an explicit expression for 
the eigenvalues in 0, it is sufficient to calculate the matrix 
elements of the operator (E — fT 22 ) _1 between the states 
\(j>kd) corresponding to a singlet formed between Fermi 
sea holes of momentum k and an electron at the impurity. 
These states are accessed by only one application of the 
different k contributions of the Hamiltonian H 12 to ) in 
the 5i subspace. 

Furthermore, as we create the space S 2 by successive 
applications of the Hamiltonian to the ground state of the 
Fermi sea, since the Hamiltonian commutes with the to¬ 
tal spin operator, the space so created contains only states 
that are singlets. Although this procedure is not capable 
of reaching the complete Hilbert space of the system, the 
ground state, our main object of analysis, is contained in it 
as it is a singlet. 

Up to this point, the procedure has been completely 
general and equal to the one presented in previous work 
done in the infinite U case E2- Now, in order to gener¬ 
alize the treatment we consider arbitrary values of U, and 
hence processes where the impurity is doubly occupied. In 


9kk (E) 


1 

E + efc — Ed — F 0 (E + efc) — F2(E + £fe) 


( 6 ) 

where the functions Fq and F 2 satisfy a set of three self 
consistent equations, exact up to terms proportional to V 2 . 


Fo(E) = Ey E-e K -F 1 (E-e K ) 

Fl(E) = 2-^r J2k E+e k -E d -F 0 (E+e k )-F 2 (E+e k ) ' 

F 2(E) =^wYjk E +e k -2E d -U-F k (E+t k ) ' ^ 

The system of equations incorporates an extra function, 
F 2 , not present in the infinite U case, related to processes 
in which the impurity is doubly occupied. 

The ground state energy, according to equations 0 and 
0 and 0, can be self consistently expressed by. 


E = F 1 (E) (8) 

The self consistent condition represents a non-perturbative 
treatment to obtain the ground state energy of the many- 
body Anderson Hamiltonian. 

Going beyond the diagonal term g kk , we incorporate 
the non diagonal elements g k k' , to the calculation of the 
ground state energy. Instead of incorporating this crossing 
terms self consistently into equation [7] which is a difficult 
numerical task, we express this non-diagonal contribution 
to the energy in terms of the auxiliary functions Fq, F\ and 
F 2 calculated as we proposed above. 

Qi(E)= 2 E £ g kk ,(E), (9) 

k,k'^k 

where the non diagonal element g k k J is 

!V\ V 2 1 

gkk'(E) = — E + _ Ed _ p ^ E + ek) _ p ^ E + x 

1 

E + efc + ek' — 2Ed — U — F\ ( E + — ek 1 ) 

1 

E + e k i - E d - F 0 (E + ek>) - F 2 (E + e k <)' 
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which gives in total, taking into account the V 2 in Eq. 0. 
a contribution proportional to V 4 to the energy. 

When Fq, F\ and F 2 are calculated through an iterative 
numerical convergent process imposed by the fulfillment 
of Eq. 0. it is possible to obtain the ground state energy 
of the system by the solution of Eq. 0 that can be written 
as, 

E = F 1 (E) + Q 1 (E) (10) 

In what follows we discuss the impurity properties de¬ 
rived by the numerical solution of the Eq. (ITOt . 

3 Numerical Results To solve the self-consistent 
equations, we consider the hybridization function to be a 
step function A(e) = A6(D — |e|) where A = ttV 2 /2D, 
and 2D is the bandwidth. From now on, we set the hy¬ 
bridization A = 1 as the unit of energy. 

Ground state energy and impurity valence 

As we state in the introduction, it is important for many 
applications to get a simple approximation that can be 
used, when the Coulomb repulsion is large enough but still 
finite, for an examination of the magnitudes derived from 
the ground state energy, such as the valence of the impurity. 

On the other hand, when the Coulomb repulsion is not 
too large, that is, of the order of the hybridization energy A, 
perturbation theory in the parameter U provides accurate 
results. For instance, renormalized-PT has been success¬ 
fully used to analyze the temperature and voltage univer¬ 
sal dependence of the conductance 0241125 1 . Specifically, 
in Fig.Q]we show the total occupancy of the impurity ob¬ 
tained from PT improved in such a way that it reproduces 
the atomic limit and satisfies the Friedel sum rule for all 
occupancies H22H23I . and the exact one calculated from the 
Bethe ansatz (see appendix) for several values of U as a 
function of the impurity energy Ed- 

As it can be seen from Fig. 0, the impurity occupa¬ 
tion obtained from PT only presents small deviations as 
compared with the exact ones (less than 1%) when the 
Coulomb repulsion is set to be U = 6A Unfortunately, 
the deviations from the exact values quickly grow in the 
strong interacting limit, when the Coulomb repulsion is 
several orders of magnitude larger than A. Also our re¬ 
sults (not shown) indicate that the magnetic susceptibility 
within the modified PT approach underestimates the mag¬ 
netic susceptibility by more than 6% for U > 3 A in the 
symmetric case Ed = —U /2. 

This situation emphasizes the convenience of having an 
approach capable of studying the strong interacting limit, 
U » A since it represents a realistic case of most semi¬ 
conductor and molecular quantum dots. The SCH expan¬ 
sion we propose to calculate the ground state energy, im¬ 
purity occupation, magnetization and susceptibility satis¬ 
fies this condition. In what follows we present our main 
results. 



Figure 1 (Color online) Total occupation of the impurity 
as a function of the energy level Ed- Solid lines stand for 
the Bethe ansatz results while symbols indicate the results 
from PT. 



Figure 2 (Color online) Total impurity occupancy as a 
function of the energy level Ed in the strong interacting 
limit for several values of the Coulomb repulsion U. Solid 
lines stand for the Bethe ansatz results. 


In Fig. 0 we show the impurity occupation as a func¬ 
tion of the local level Ed for several values of U calculated 
from the SCH expansion including the 4th order processes. 
The solid lines indicate the results obtained from the Bethe 
ansatz technique. We include an appendix in which we ex¬ 
plicitly show the steps we follow to compute the integral 
expressions that determine the valence within the Bethe 
ansatz solution for arbitrary values of the model parame- 
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Figure 3 (Color online) Universal occupation number rid 
as a function of the scaled energy nE^/2A for the same 
parameters that in Fig. ©. 

ters. Here, we have changed the unit of energy A = 1 by 
setting A = 0.1; 0.2 and 0.3 while we chose U = 8. 

In Fig. © we plot the SCH results for the impurity 
occupancy shown in Fig. © as a function of the renor¬ 
malized energy E% = Ed + —1 n(^^) which explicitly 
includes the Haldane shift together with the corre¬ 
sponding renormalized Bethe ansatz curve. As expected 
from the universality of the model, we notice here that 
all the curves collapse into an universal one in the Kondo 
regime —E^^$>A and that at E* t =0 the valence is rid = 
0.58 in agreement with Fig. 2 of Ref. ©. 

Magnetic properties 

In what follows, we add a magnetic field B to the sys¬ 
tem and analyze the magnetization of the impurity, Mi mp . 
The inclusion of the Zeeman interaction within the SCH 
procedure is straightforward and it is described in Ref. 
If20ll . According to the usual definition, I©, to obtain Mi mp 
we subtract from the total magnetization of the system the 
magnetization of the isolated Fermi sea. Therefore, the im¬ 
purity contribution to the magnetization, Mi mp , is given 

by, 

M imp (B) = —^(E(B)-e T (B)). (11) 

In Fig. ©, we show the impurity contribution to the 
magnetization in units of <jfin, as a function of B. We 
chose Ed = —6 to compare our results with the exact 
Bethe ansatz results for the Kondo model |[28ll . The se¬ 
lected values of the energy level Ed = —6 and U = 12 
in units of A, corresponds to the symmetric case, in which 
the charge fluctuations are frozen and the Anderson model 
maps onto the Kondo one. 



Figure 4 (Color online) Comparison of the results ob¬ 
tained by the SCH expansion for the symmetric case and 
the Bethe ansatz solution of the Kondo model, for the im¬ 
purity contribution to the magnetization in units of 

g^B as a function of B/T\. The energy impurity level is 
taken to be Ed = —6. The scale T\ is defined in the main 
text. 


Following the notation introduced by N. Andrei in Ref. 
f28ll . we define the energy scale T), that characterizes the 
strong coupling regime (here it means Kondo regime) 126S . 
from the susceptibility at zero temperature, 

Ximp (B = 0) = (12) 

J i 

It can be noticed from Fig. © that the agreement be¬ 
tween SCH expansion and BA is very good when Mi mp 
is calculated in units of B/Ti , for which Mi mp is an uni¬ 
versal curve. For small values of B as compared with 7 j, 
we obtain a linear dependence of the magnetization with 
respect to the applied field. This agrees with the expected 
behavior of the Kondo problem, which is indicative of a 
screened impurity spin. On the other hand, for larger fields, 
we obtain the characteristic slow asymptotic approach to 
the value g/rs/2 of saturation, which also agrees with the 
well known logarithmic corrections lf28l [3l |29l . 

From the definition of the low energy scale T), it is 
natural |[26l to r elate it with the usual Kondo scale, Tk 
via T\ = \Z8ne Tk ■ Unfortunately, the Kondo scale that 
we obtain within the SCH procedure is found to be several 
orders of magnitude smaller than the expected one from 

the well known Haldane scale, Tk = \J^-exp(—^) 
m, for the symmetric case. For the parameters used in 
Fig. [U the SCH expansion predicts a Kondo scale of the 
order of 10 -5 while the corresponding one to the Haldane 
expression is of the order of 10 -2 in units of A. 

This is not a surprising result. The same feature appears 
when the SCH expansion in its lowest version, NCA, is ap- 
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Figure 5 Impurity contribution to magnetic susceptibility 
X* = x/ (i^Bg) 2 /A)as a function of the renormalized en¬ 
ergy Ej. The solid line indicates the calculation by using 
the Bethe ansatz technique in the limit of both, D, U —> oo. 
The symbols represent the SCH solution for two selected 
values of the bandwidth D. 


plied using finite values of the Coulomb interaction 111211131 
HU. As we mentioned in the introduction, this pathology 
was partially overcome, adding in a self consistent way, 
all the dressed diagrams that involve one crossing conduc¬ 
tion electrons, which corresponds to the next leading or¬ 
der called one-crossing approximation OCA. While in our 
present treatment of the problem we are already taking into 
account 4th order processes, we are not including these 
contributions in a self consistent way, as we mentioned in 
section [2] and can be seen from Eq. l[7J). While our treat¬ 
ment is completely reliable to obtain the static properties of 
the Anderson Hamiltonian in the strong coupling limit, it 
fails to give correctly the Kondo temperature when U/A is 
not large enough, situation in which an equivalent scheme 
to OCA approximation is required. This full treatment is 
beyond the scope of the present study. 

As a final analysis, we include in Fig. El the suscepti¬ 
bility comparing our results for Ximp with those obtained 
from a Bethe ansatz calculation l30l in the case of infinite 
Coulomb repulsion. In this case, we show results for both 
D = 10 and D —> oo properly scaled by a Haldane shift. 
The agreement with the Bethe ansatz solution of Ximp is 
very good in both cases, capturing the universality of the 
model. 


previous treatment for infinite U. We show that the static 
properties of the model can be well calculated within the 
SCH expansions for zero temperature and finite magnetic 
fields. We successfully match our results with those com¬ 
ing from the Bethe ansatz solution of the model. 

This simple approach allows the study of systems that 
can be described by the Anderson model, such as semi¬ 
conductor EQ or molecular quantum dots |[32l where 
the Coulomb repulsion is expected to be large but cannot 
be considered infinite. Experimentally relevant quantities, 
such as the equilibrium conductance of the impurity, which 
is directly related to the occupancy (if the system is in a 
Fermi liquid regime) can be obtained. Furthermore, the 
semi analytic nature of the solution permits a direct and 
fast calculation of the properties on a single workstation. 

A Bethe ansatz equations of the impurity valence 
at zero magnetic field To compute the valence of the 
impurity at zero magnetic field and for arbitrary values of 
the parameters U, e d , A we start with the expressions (5.8) 
and (5.9) given by Wiegmann and Tsvelick in Rcfs|j2][3|. 


U/2 + e d _ i f + °° dco e ~^/ 2 ~ iuiQ 1 
(2 UA) 1 / 2 v^Z-oo uJ + ir] G(“)(w) (-iu + r ]) 1 / 2 

(13) 

i r + °° du e-H/ 2 

n d = 1 - 


v^W-oo w + M7G(-)(w) 


r»+oo 




(14) 


where the functions G 1 ' \uj), g(n) and Z\(«) are given by 


G*- 1 (w) = V2tt 



ff(«) 

A(k) 


(K-e d - U/2) 2 
2 UA 
A 1 

7T (ft - Crf) 2 + Z\ 2 ’ 


with ry —> 0 + . The equation in (fl3l > fixes the value of Q that 
enters in the explicit calculation of the impurity occupation 
in (fril l. 

The highly divergent nature of the integrand in Eq. 
© demands particular care in handling the limit rj —> 
0+. In fact, using the identity —— = —)=— ' - —> 

D J ^y—iuj+T) sj—i ^/ui+ip 

-— (1 + i sgn{u> )) and taking explicitly the limit g —> 
v 2 M 

0 + the RHS of O- which we call Z(Q), becomes 


4 Summary and Conclusions In this work we ap¬ 
plied the self consistent hybridization expansions including 
4th order processes to compute the ground state energy, oc¬ 
cupancy and magnetic properties of the Anderson impurity 
model in the case of finite Coulomb repulsion, extending a 


Z(Q) = 


J-f 
y/2^ J- 

..J-I 

V2n J- 


+ oo 


duj 


o-H/ 2 -iuQ 


to + ig G( _ )(w) (—iw + ry) 1 / 2 

+ °° duJ n / • At 

( 1 /CJ — I 7 vd{i 0 )) 


\/2R 


x (1 + i sgn(u})) H(u, Q), 
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with H (at, Q) 


e -M/2-iuQ ! 


72 


M74 




and being the limit 77(0, Q) = 1 /%/2, the argument has a non 
integrable singularity at a; = 0. In order to avoid this singular¬ 
ity and make the integral tractable numerically, we subtract the 
singular part of Z(Q) before taking the limit 77 —> 0 + . 

Let us write the integral Z(Q) in the following form 


Z(Q) 


i 

\J—2m 



dot e —M/2 -iu,Q 

(a j + i V ) 3 / 2 G(-)(at) 


(15) 


Using the fact that 




ic 

\/—2m 



duj 

(cj + irf) 3 / 2 


= 0 


we set c = 1/G*' ^(0) = l/\/2 and re-write Z(Q) 
as follows 


(16) 
Z(Q) - 


Z(Q) = 


= t 

■2iir J _ 


’ / 


, da; ,, /2 { H(u,Q)~ 4= } 

(w + IT]) 3 / 2 1 ' ^2 


20 F y_ 


— 00 


CJ 


vr 


72 


(17) 


In this way, the argument within the integral becomes regular. 
We also map the infinite range of integration to a finite one 
by using 


/ +00 r+oo 

duj /(or) = duj {/(w) + /(-oj)} 

-00 «/ 0 


= f X 

Jo 


(1 — x) 3 (1 — x)‘ 




The change of variable was chosen so that the argument of 
the integral is smooth for x —> 1 (u> —> 00 ). By means of the 
above mentioned transformations, the integral becomes suitable 
for numerical evaluation. 

The integral in GS is simpler and it can be evaluated in a 
more direct way. Calling K(uj) = e »“( 9 ( K )-< 5 )/\( / j( 

becomes 


rid = 1 — 


i r+°° 

y/2n J — 00 


duj 


o l w l /2 


V2n J-ao w + ir\ G( ((w) 


K( u 


= 1 - 


= 1 - 


J-f 

72 tt y_ 


duj ( 1 /a; — iir6(uj)) 


e -l“l/2 

G(-)(a;) 




72: 


1 


_ 1 

_ 2 + V^TT .L 


b 7 
/_ 


+00 g-l^l/2 j 

-00 


G( _ )(ai) 


72 J 


+oo g-l-l /2 K{uj) 

a; G^ ■'(ai) 


(18) 


Finally, we compute the principal part of the integral with the 
help of the following transformations 


/ + 00 c+oo 

dx> f(ui) = du {f(uo) + /(-oj)} 

-00 J 0 

R dI (TR {/( T7 )+/( 77- 
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